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CHAPTER  I 


INTRODUCTION 

The  value  of  the  Impulse  response  waveform  as  a  means  of 
characterizing  a  radar  target  was  suggested  by  Kennaugh  and  Cosgrlff  [1] 
In  1958  and  refined  by  Kennaugh  [2],  and  Kennaugh  and  Moffatt  [3],  [4]. 
The  Impulse  response  Is  that  function  of  time  which  describes  the  field 
scattered  by  a  radar  target  Illuminated  by  an  Impulsive  plane  wave.  The 
impulse  response,  Fj { t) ,  and  the  normalized  transfer  function,  F(jw),  of 
the  scatterer  are  a  Fourier  transform  pair.  The  Ideal  Impulse  response 
can  only  be  approximated,  whether  by  theory  or  experiment,  and  while 
such  approximations  show  varying  degrees  of  success,  the  term  Impulse 
response  will  be  used  herein  to  refer  to  the  entire  range  of 
approximations. 

The  value  of  the  Impulse  response  lies  mainly  In  Its  Intimate 
relationship  to  the  geometry  of  the  target.  The  times  of  arrival  of  the 
contributions  from  Individual  scattering  mechanisms  on  the  target  can  be 
easily  predicted  on  the  time  axis  of  the  impulse  response  waveform. 

This  relationship  provides  an  Intuitive  "feel"  for  the  scattering 
properties  of  each  target  geometry.  The  Insight  It  affords  Is 
especially  valuable  when  used  in  conjunction  with  high  frequency 
scattering  theories,  such  as  the  Uniform  Theory  of  Diffraction  (UTD), 
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which  localize  and  superpose  the  contributions  from  Individual 
scattering  mechanisms  on  the  target. 

The  early  studies  mentioned  above  used  a  combination  of  time  and 
frequency  domain  techniques  and  low  and  high  frequency  solutions  to  the 
scattered  field  to  assemble  Impulse  response  approximations  for  several 
simple  geometries.  The  Rayleigh  law  of  scattering  provided  a  set  of 
moment  conditions  on  the  Impulse  response,  while  physical  optics 
predicted  the  high  frequency  behavior.  The  approximations  were  greatly 
Improved  by  using  the  high  frequency  predictions  of  the  Geometrical 
Theory  of  Diffraction.  Various  estimates  of  the  Impulse  response 
waveform  for  several  simple  targets  are  derived  In  [4]  and  dlcussed  In 
more  detail  In  Chapter  IV. 

These  studies  sparked  an  Interest  In  coherent  frequency  domain 
backscatterlng  measurements,  since  |F(jw)l2  Is  related  to  the  radar 
cross  section  of  a  scatterer  by  a  multiplicative  constant.  A  measured 
Impulse  response  waveform  may  be  obtained  as  the  Inverse  Fourier 
transform  of  F(ju)  measured  over  a  broad  band  of  frequencies.  The 
bandwidth  requirement  Is  the  primary  problem  In  the  measurement  of 
Impulse  responses.  Approximations  of  the  ramp  response,  the  second 
Integral  of  the  Impulse  response,  have  been  constructed  from 
measurements  at  just  ten  harmonically  related  frequencies  [5].  While 
the  ramp  response  is  valuable  In  such  radar  signature  areas  as  target 
Imaging  [63,  It  will  not  yield  the  impulse  response.  A  bandwidth  of 
approximately  100:1  Is  necessary  to  obtain  a  good  Impulse  response 
[63. 

2 


1 


The  present  study  extends  the  measured  transfer  function  to  a 
practically  continuous  spectrum  covering  a  96:1  bandwidth.  Using  target 
dimension  scaling  on  a  broadband,  stepped  frequency  measurement  system. 
It  was  possible  to  measure  F(jui)  at  frequencies  beginning  In  the 
Rayleigh  region,  spanning  the  resonance  region,  and  extending  Into  the 
optics  region.  This  Is  the  first  attempt  to  extend  the  bandwidth  of  the 
measured  F(jw)  to  that  sufficient  for  the  construction  of  an  accurate 
approximation  to  the  measured  Impulse  response  waveform. 

One  of  the  objectives  of  this  study  has  been  to  assemble  "complete" 
complex  scattering  cross  section  data  for  several  canonical  radar 
targets.  The  data  are  complete  In  the  following  respects: 

(1)  Both  amplitude  and  phase  of  the  backscattered  signal  were 
measured. 

(2)  The  data  cover  a  96:1  bandwidth. 

(3)  The  backscattered  signal  was  measured  for  three  polarizations 
of  the  transmit  and  receive  radar:  vertical -vertical , 
horizontal-horizontal ,  and  vertical-horizontal . 

(4)  The  measurements  were  made  at  look  angles  of  every  15  degrees 
from  0  to  180  degrees. 

(5)  The  vertical -horizontal  (cross)  polarized  measurements  were 
made  at  aspects  out  of  the  principal  planes  of  the  targets. 

The  data  bank  assembled  from  these  measurements  is  of  Interest  to 
those  doing  Inverse  scattering  research,  to  scattering  analysts,  and  to 


those  interested  in  the  use  or  control  of  a  radar  target's  scattering 
properties. 

Chapter  II  describes  the  frequency  domain  range,  equipment,  and 
target  geometry;  Chapter  III  discusses  the  data  processing  utilized; 
Chapter  IV  presents  measured  results  and  compares  those  results  to 
theoretical  solutions;  Chapter  V  discusses  the  cross  polarized 
measurements.  Finally,  the  relation  between  the  impulse  response 
waveform  and  the  high  frequency  scattering  analysis  of  the  Uniform 
Theory  of  Diffraction  is  investigated,  and  the  potential  for 
complementary  interaction  between  the  two  is  demonstrated. 
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CHAPTER  II 


BROADBAND  BACKSCATTERING  MEASUREMENT  SYSTEM 

The  Impulse  response  may  be  measured  on  a  time  or  a  frequency 
domain  range.  A  time  domain  range  has  two  disadvantages  which  restrict 
the  ^solution  of  the  measured  Impulse  response.  The  first  Is  the 
difficulty  of  generating  an  Impulse  with  sufficient  power  and  short 
enough  duration  ( <20  psec)  to  achieve  adequate  spectral  coverage. 
Assuming  such  an  Impulse  can  be  generated,  the  second  problem  Is  the 
need  for  a  single  broadband  antenna  to  transmit  it. 

A  frequency  domain  range  has  three  advantages  over  the  time  domain 
range,  each  of  which  hinges  upon  subdividing  the  total  band.  First,  the 
required  bandwidth  need  not  be  covered  by  a  single  antenna.  Second, 
scaling  can  be  used  to  Increase  the  available  bandwidth  of  the  system 
hardware  by  the  scale  factor  of  the  models.  Third,  the  power  capability 
of  the  radar  system  Is  greater  over  a  restricted  band.  To  Implement 
this  frequency  band  subdivision  requires  a  certain  amount  of  hardware 
flexibility,  but  It  can  be  achieved  with  off-the-shelf  microwave 
components. 

The  obvious  disadvantage  Is  that  more  measurements  need  be  made  In 
the  frequency  domain,  and  both  amplitude  and  phase  need  to  be  stored  at 
each  frequency.  A  computer  controlled  RCS  measurement  system  which  can 
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speed  up  data  acquisition  makes  the  measurement  of  accurate  Impulse 
response  waveforms  feasible. 

A  broadband,  frequency  domain,  complex  radar  cross  section 
measurement  range  based  on  the  concept  of  background  nulling  has  been 
Implemented  for  this  purpose.  The  system  Is  similar  to  that  devised  by 
Weir  et  al.,  at  Stanford  [18].  By  circumventing  the  need  for  a  time 
consuming  background  nulling  at  each  frequency  the  system  allows  for  the 
rapid  collection  of  target  backscatter  over  a  wide  band  of  frequencies. 
The  subtraction  and  calibration  algorithm  will  be  explained  In  detail  In 
Chapter  III. 

The  data  collection  system  Is  shown  In  Figures  2-1  and  2-2. 

Data  acquisition  Is  controlled  by  a  Declab  PDP11/03  minicomputer.  The 
source  Is  a  Hewlett  Packard  8690B  microwave  oscillator  with 
Interchangeable  octave  band  modules  of  1-2  GHz,  2-4  GHz,  4-8  GHz,  and 
8-12  GHz.  No  external  amplification  Is  used  and  output  power  varies 
from  10-70  mW  depending  on  the  module.  Before  each  reading  Is  taken  an 
Autohet  microwave  counter  checks  the  frequency  of  the  output  and 
corrects  It  to  within  a  resolution  of  1  MHz. 

The  receiver  Is  a  Hewlett  Packard  8410A  Network  Analyzer.  The 
received  signal  Is  a  vector  measured  relative  to  a  reference  signal 
tapped  off  the  transmit  port.  The  Network  Analyzer  provides  relative 
amplitude  and  phase  (50  mV/dB  and  10  mV /degree)  as  analog  voltages  on 
separate  channels.  The  reference  signal  Is  routed  through  a  delay  line 
which  approximates  the  average  two-way  phase  path  of  the  anecholc 
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Figure  2-1.  Broadband,  frequency  domain,  complex  radar  cross 
section  measurment  range. 


chamber.  This  delay  minimizes  the  number  of  phase  branch  cuts  In  the 
received  signal  relative  to  reference. 

The  measurement  facility  uses  two  broadband  AEL  microwave  antennas 
mounted  back  to  back  at  the  focal  point  of  a  three  foot  parabolic 
reflector.  Such  an  arrangement  allows  the  measurement  of  true 
backscatter.  The  antennas  and  the  range  are  shown  In  Figure  2-3.  The 
transmit  and  receive  horns  are  Isolated  with  absorbing  material  to 
reduce  direct  coupling. 

While  It  Is  not  necessary  to  null  the  background  with  this  system. 
It  Is  desirable  to  reduce  the  strength  of  multipath  returns,  such  as  the 
signal  scattered  from  the  target  to  the  floor  and  reaching  the  receive 
antenna.  For  this  reason  the  chamber  Is  kept  as  anecholc  as  possible. 
The  background  return.  Including  antenna  coupling,  was  reduced 
approximately  60  dB  below  the  direct  coupled  signal. 

The  AEL  microwave  antennas  operate  from  1  to  12  GHz.  The  targets 
were  constructed  of  cast  aluminum  In  two  scale  sizes  related  by  an  8:1 
scale  factor.  This  combination  gives  an  effective  bandwidth  of  96:1, 
which  Is  In  the  range  required  for  accurate  Impulse  response 
construction.  The  targets  of  primary  Interest  In  this  study  are:  a  2:1 
right  circular  cylinder,  a  2:1  right  circular  cylinder  with  a 
hemispherical  cap  on  one  end,  and  a  2:1  prolate  spheroid.  Target 
geometry  and  actual  physical  dimensions  are  shown  In  Figure  2-4.  The 
frequency  ranged  over  0.9  <  kL  <  96,  where  L  is  the  maximum  target 
dimension  and  k  *  2*/x.  Spheres  ranging  In  diameter  from  1.25  to 
6.0  Inches  were  measured  over  the  same  frequencies. 


Figure  2-3.  Anecholc  chamber  and  back  to  back  microwave 
antennas  set  up  to  measure  true  backscatter 


Figure  2-4.  Target  set: 

a)  2:1  right  circular  cylinder  of  radius  a. 

b)  2:1  right  circular  cylinder  of  radius  a  with 
hemispherical  cap. 

c)  2:1  prolate  spheroid. 

d)  calibration  targets  and  scale  targets. 
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The  targets  were  mounted  on  a  ten  foot  styrofoam  pedestal  which  Is 
rotated  under  computer  control.  Measurements  were  made  at  increments  of 
fifteen  degrees  in  both  the  E  plane  and  the  H  plane.  For  the  cross 
polarized  measurements  the  target  set  included  discs,  flat  plates,  and 
strips  at  look  angles  out  of  the  planes  of  symmetry.  The  cross  polarized 
measurements  are  discussed  separately  in  Chapter  VI. 


CHAPTER  III 


DATA  PROCESSING 

The  discussion  of  data  processing  covers  the  following  areas: 

1)  calibration  of  the  measured  data,  2)  sampling  criteria,  3)  smoothing, 
4)  windowing,  5)  resolution  In  the  time  domain,  6)  possible  sources  of 
error,  and  7)  error  estimation. 

The  calibration  algorithm  is  given  In  Equation  (3.1).  The 
broadband  measurement  range  uses  background  subtraction  rather  than 
background  nulling.  The  amplitude  and  phase  of  the  backscattered  signal 
Is  recorded  at  discrete  frequencies  for  each  of  three  cases:  1)  target 
plus  background,  2)  background  only  with  no  target,  3)  calibration 
target  plus  background.  The  calibration  target  Is  one  with  a  known 
exact  complex  backscatter  over  the  frequencies  of  Interest.  For  the 
co-pol arl zed  measurements,  spheres  of  various  sizes  are  used  as 
calibration  targets. 

The  measurments  of  target  plus  background  (TB)  are  stepped  through 
a  set  of  pre-determl ned  frequencies.  Normally  201  data  points  are 
recorded  over  an  octave  band.  The  measurements  of  background  with  no 
target  (B)  are  then  stepped  through  the  same  frequencies.  The  process 
Is  repeated  for  the  calibration  sphere  (CB).  At  each  frequency  these 
three  sets  of  measured  data  are  then  combined  with  the  exact  sphere 
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spectrum  (E)  In  the  following  algorithm: 


Tc  =  (TB-B)  .(E)  _ 
(CB-B) 


(3.1) 


Since  the  factor  (E )/( CB-B)  Is  a  complex  calibration  constant  which 
calibrates  out  the  system  variables,  the  result,  Tq,  Is  the  complex 
radar  cross  section  of  the  target  alone.  The  problem  of  noise  Is 
considered  In  subsequent  paragraphs. 

When  spectral  data  are  sampled  and  passed  through  a  Fast  Fourier 
Transform  (FFT),  the  resulting  time  waveform  Is  periodic.  Overlap,  or 
aliasing.  In  the  time  waveform  Is  a  possible  source  of  error  which  Is 
avoided  by  proper  selection  of  the  frequency  sampling  Interval,  Am.  If 
the  time  domain  transient  response  Is  zero  for  lt|  >  T,  then  the  time 
waveform  Is  accurately  reconstructed  by  frequency  samples  at  Intervals 
of  am  <  */T.  It  Is  safe  to  assume  that  the  transient  response  dies  out 
after  T  »  6t,  where  t  ■  L/C,  L  Is  the  maximum  target  dimension,  and  C  Is 
the  speed  of  light.  Thus  the  frequency  spectrum  need  only  be  sampled  at 
Intervals  of  Am  <  V6T,  or  AkL  <  V6,  where  Ic  «  Zk/X. 

The  measured  spectral  data  are  considerably  oversampled,  with  the 
sampling  rate  ranging  from  five  times  to  one  hundred  times  the  minimum 
rate.  A  total  of  1401  unequally  spaced  frequency  samples  are  distilled 
prior  to  the  FFT  to  kL  samples  at  Intervals  of  AkL  ■  1/10,  which  Is  5.24 
times  the  minimum  rate.  This  sampling  Interval  fixes  the  total  time 
frame  at  63 t. 


Oversampling  justifies  smoothing  to  Improve  the  accuracy  of  the 
measured  spectral  data.  At  each  sample  point,  a  smoothed  phasor  value 
Is  found  by  weighting  the  real  and  Imaginary  parts  of  the  (N-l)/2 
adjacent  samples  on  either  side  of  It  with  a  shifted  cosine.  The 
weighting  function  Is 

w(n)  =  0.5[1.0  +  cos(2nif/M)],  n=  -(N-l)/2,  ...»  0,  ...(N-l)/2  , 

(3.2) 

with  the  sum  of  the  weights  equal  to  N. 

The  number,  N,  of  samples  averaged  varies  from  5  to  25  depending  on 
frequency.  Moving  the  weighting  function  across  the  frequency  samples 
Is  a  convolution  In  the  frequency  domain.  This  Is  equivalent  to 
multiplication  In  the  time  domain  by  the  Inverse  transform  of  w(n), 
which  Is  shown  In  Figure  3-1  for  N*15.  The  target  information  lying  In 
a  6t  time  frame  centered  at  t=0  Is  passed  unattenuated,  while  the  range 
clutter  lying  outside  ±6t  Is  effectively  removed.  Smoothing  thus  serves 
to  Improve  the  accuracy  of  the  samples  through  averaging,  as  well  as  to 
time  gate  some  undesired  signals. 

/w  (V 

In  Equation  (3.1)  the  difference  phasors,  (TB-B)  and  (CB-B),  both 
contain  noise.  It  has  been  found  that  smoothing  both  difference  phasors 
before  division  Inhibits  the  multiplication  of  noise  Inherent  In  the 
division.  This  operation  requires  that  the  range  phase,  e"J2kR,  be 
factored  out  of  both  numerator  and  denominator,  or  equivalently,  that 
the  difference  vectors  be  centered  In  the  time  window  shown  In  Figure 
3-1.  Some  pre-smoothing  of  the  phase  has  also  been  found  to  Improve  the 


Figure  3-1.  Inverse  Fourier  transform  of  the  frequency 
domain  weighting  function,  1.  +  1,  cos(2n*/N), 

for  N  ■  15.  Unit  of  t*L/C  where  L  ■  maximum 
target  dimensions,  and  C  ■  speed  of  light. 


results,  especially  in  those  parts  of  the  spectrum  where  phase  branch 
cuts  occur. 

Windowing  of  the  frequency  domain  data  is  used  to  reduce  the  Gibb's 
phenomenon  caused  by  truncation  of  the  measured  spectrum.  Hanning, 
Hamming,  Gaussian,  and  Tukey  windows  have  been  implemented  for  this 
purpose.  The  Hanning  window  has  been  found  to  give  the  best  overall 
results.  The  weighting  function  for  the  Hanning  window  is  the  same  as 
that  given  In  Equation  (3.2)  except  that  the  weighting  is  applied  across 
the  entire  spectrum.  Such  multiplication  In  the  frequency  domain  is 
equivalent  to  convolution  in  the  time  domain.  The  Inverse  transforms  of 
the  various  windows  are  given  in  [19]  along  with  their  performance 
characteristics.  In  particular,  the  highest  sidelobes  of  the  transform 
of  the  Hanning  window  are  -32  dB  as  compared  to  -13  dB  for  the  sin(x)/x 
transform  of  the  rectangular  window.  The  improvement  in  si  delobe  level 
is  in  exchange  for  a  50  percent  Increase  in  the  equivalent  pulse  width 
from  0.066  t  for  the  rectangular  window  to  0.10  t  for  the  Hanning 
window.  The  Inverse  transform  of  the  Hanning  window  Is  shown  in  Figure 
3-2  for  kL  <  96.0.  On  the  linear  amplitude  scale  the  sidelobes  are 
barely  discernable. 

The  equivalent  pulse  width  discussed  above  Is  a  measure  of  the 
resolution  of  the  time  domain  waveforms.  The  resolution  is  fixed  by  the 
maximum  frequency  of  the  samples.  For  the  15  inch  target  we  will 
examine  In  Chapter  IV,  t  =  1.27  nsec,  and  for  a  maximum  sampling 
frequency  of  12  GHz  the  equivalent  pulse  width  Is  0.0833  nsec, 
increasing  to  0.127  nsec  when  the  Hanning  window  is  applied.  The  time 
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scale  In  Figure  3-2  may  be  used  as  a  reference  for  the  resolution  In  the 
waveforms  to  be  presented  In  Chapter  IV. 

Before  moving  on  to  a  discussion  of  noise  one  last  point  Is  worth 
noting.  Since  the  measurements  extend  well  Into  the  Rayleigh  region  It 
Is  possible  to  apply  a  least  squares  fit  to  the  first  200  samples,  from 
which  a  Rayleigh  coefficient  may  be  extrapolated.  This  allows  the 
spectrum  to  be  extended  down  to  the  first  harmonic  at  kL  =  0.1,  thereby 
fitting  the  lower  harmonics  and  eliminating  the  discontinuity  at  kL  = 
1.0. 

Besides  those  already  mentioned,  other  possible  sources  of 
measurement  error  are  system  noise,  multi -path  scattering,  and  target 
positioning  error.  System  noise  Includes  amplifier  noise  In  the 
receiver  front  end  and  In  the  A-D  converters  and  frequency  drift.  A 
figure  of  merit  for  system  noise  is  the  noise  floor.  The  noise  floor  Is 
measured  by  treating  a  background  measurement  as  a  target  in  Equation 
(3.1)  and  subtracting  another  background  measurement,  taken  during  the 
same  data  run.  The  result  Is  essentially  a  repeatability  test.  The  two 
background  measurements  are  taken  at  the  beginning  and  end  of  a  data 
session  so  as  to  Include  the  effects  of  drift.  A  typical  noise  floor  of 
between  0  and  5  dB>cm2(-40  to  -35  dBsm)  is  shown  In  Figure  3-3.  The 
resulting  phase  should  be  random. 

Error  due  to  multi -path  scattering  from  the  target  to  some  part  of 
the  chamber,  and  finally  to  the  receive  antenna  Is  reduced  by  making  the 
chamber  as  anecholc  as  possible.  The  background  scattering  was  reduced 
approximately  50-60  dB  with  absorbing  material.  In  particular,  the 
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direct  coupling  between  transmit  apd  receive  horns  was  reduced  by 
wrapping  them  with  absorber,  as  can  be  seen  In  Figure  2-3. 

Target  positioning  error  may  corrupt  the  measured  phase  data.  It 
has  been  found  that  the  positioning  Is  accurate  for  frequencies  up  to  8 
GHz,  and  adjustments  need  be  made  only  In  the  8-12  GHz  band.  The 
correcting  procedure  Is  to  Inverse  transform  the  measured  8-12  GHz 
spectra  separately  before  combining  with  the  lower  frequency  data.  If 
the  scattering  from  the  specular  point  Is  used  as  a  reference  on  the 
time  axis  the  target  can  be  "moved"  to  proper  alignment  position  by  a 
simple  time  shift.  This  Is  equivalent  to  matching  the  slope  of  the 
phase  at  the  8  GHz  junction. 

Finally,  an  estimate  of  the  total  measurement  error  has  been 
performed  to  place  reliability  bounds  on  the  data.  Applying  a  single, 
one-dimensional  error  figure  to  broadband,  complex  data  is  not  valid 
since  the  measurement  error  varies  with  target  cross  section,  frequency, 
and  coherency  of  the  amplitude  and  phase  data. 

In  an  early  attempt  at  placing  error  bounds  on  the  data  a  large 
number  of  samples  were  taken  at  a  single  frequency  for  each  of  the  terms 
In  Equation  (3.1),  and  run  through  the  calibration  algorithm.  Standard 
deviations  about  the  means  were  calculated  for  amplitude  and  phase 
separately.  The  procedure  was  repeated  for  targets  with  cross  sections 
ranging  from  40  dB>cm2  down  to  10  dB>cm2,  and  at  several  frequencies. 

The  resulting  error  estimates  ranged  from  ±0.5  dB  and  ±4.0°  for  target 
cross  sections  above  +25  dB>cm2  to  ±1.7  dB  and  ±10°  for  target  cross 
sections  of  +10  dB>cm2.  It  should  be  noted  that  these  figures  are  for 
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unsmoothed  data  and  thus  do  not  reflect  the  processing  discussed 
earlier.  These  early  estimates  are  quoted  In  [20]. 

A  more  accurate  error  estimate  Is  obtained  In  the  time  domain. 
Besides  treating  amplitude  and  phase  coherently,  a  time  domain  analysis 
Includes  the  effect  of  smoothing  and  all  possible  sources  of  error.  In 
this  analysis  the  measured  data  for  a  sphere  are  compared  to  exact 
sphere  data.  There  are  several  reasons  for  using  the  sphere  error  as 
representative  of  all  measurement  error.  The  sphere  is  used  to 
calibrate  all  the  co-polarlzed  measurements;  the  cross  section  of  the 
sphere  Is  lower  than  that  of  the  larger  targets  of  Interest;  and  the 
sphere  has  an  exact  spectrum  from  which  to  calculate  measurement 
deviation. 

Figure  3-4  shows  the  exact  and  measured  spectra  for 
the  sphere.  The  measured  spectrum  Is  offset  +10  dB.  A  Gaussian  taper 
has  been  applied  to  the  highest  harmonics  to  eliminate  the  sharp 
truncation  at  12  GHz,  but  no  windowing  has  been  done  across  the  midband. 
The  corresponding  phase  Information  Is  shown  In  Figure  3-5.  Both 
spectra  are  Inverse  transformed  and  the  difference  taken  between  the 
time  waveforms.  This  difference  Is  then  transformed  to  give  the  error 
spectrum  In  Figure  3-4.  The  vertical  scale  Is  in  dB  below  the  physical 
optics  cross  section  for  the  sphere.  The  -20  dB  level  corresponds  to  1 
percent  error  or  ±0.04  dB.  It  can  be  seen  that  throughout  most  of  the 
band  the  measurement  error  falls  below  the  1  percent  level.  The  worst 
case  error  Is  at  -12  dB  corresponding  to  a  5  percent  error  or  ±0.25  dB. 
This  occurs  at  frequencies  below  1  GHz  and  for  one  peak  at  3  GHz. 
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Figure  3-4.  Upper  curve:  measured  sphere  spectrum  offset  +10  dB. 
Middle  curve:  exact  sphere  spectrum. 

Lower  curve:  spectrum  of  difference  between  measured 
and  exact. 

Vertical  scale:  radar  cross  section  In  dB  >  wa*  for 
a  *  sphere  radius  *  7.62  cm. 

Horizontal  scale:  frequency  In  GHz. 
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Figure  3-5.  Measured  versus  exact  phase  of  sphere  with  radl 
a  ■  7.62  cm. 
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In  the  low  frequency  region  the  error  Is  due  to  the  generally  lower 
cross  section  of  the  targets  combined  with  the  difficulty  of  measuring 
phase  as  It  hovers  around  180°.  At  3  GHz  the  amplitude  spectra  compare 
very  well,  but  Figure  3-5  shows  a  phase  error  near  a  branch  cut.  This 
Illustrates  that  the  time  domain  error  analysis  Includes  the  effect  of 
phase  coherently. 

To  demonstrate  that  the  Improved  error  performance  shown  In  Figure 
3-4  Is  due  to  smoothing  we  may  compare  It  with  the  error  estimates  for 
unsmoothed  data  given  earlier.  One  would  expect  the  standard  deviation 
for  the  unsmoothed  data  to  be  reduced  by  a  factor  where  N  Is  the 

number  of  samples  Included  In  the  smoothing.  A  standard  deviation  of 
0.5  dB  In  the  unsmoothed  data  corresponds  to  an  error  of  22  sq.  cm.  for 
a  six  Inch  diameter  sphere  with  physical  optics  cross  section  of  182  sq. 
cm.  If  a  21  point  smoothing  Is  applied,  this  error  Is  reduced  to  4.85 
sq.  cm.,  which  Is  16  dB  down  from  the  physical  optics  cross  section. 

This  figure  Is  closer  to  the  error  estimate  In  Figure  3-4. 

In  summation,  the  measured  data  Is  accurate  to  within  1  percent  or 
better  over  most  of  the  spectrum.  This  figure  can  be  applied  to  phase 
and  amplitude  together,  and  Includes  all  possible  sources  of  error. 
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CHAPTER  IV 


MEASURED  IMPULSE  RESPONSES 


The  Impulse  response  of  a  scattering  object  Is  a  time  domain 
concept  derived  from  linear  systems  analysis  [1].  With  the  input 
defined  to  be  the  linearly  polarized  field  Incident  upon  the  scatterer 
the  output  Is  taken  to  be  a  selected  component  of  the  scattered  field  in 
a  particular  direction.  The  derivation  of  the  Impulse  response  concept 
applied  to  E.M.  scattering  will  be  briefly  summarized  here  as  will  the 
physical  optics  approximation  to  the  impulse  response  of  a  few  simple 
shapes. 

The  radar  cross  section,  a,  of  a  scattering  object  Is  a  frequency 
domain  concept  defined  as 
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ES(Ju) 

E*'(Ju) 


(4.1) 


where  Es(ju>)  Is  a  selected  polarization  component  of  the  scattered 
electric  field  at  a  distance  r  from  the  object  and  E^(ju)  Is  the 
linearly  polarized  electric  field  Incident  upon  the  object.  The 
geometry  of  the  scattering  problem  Is  sketched  In  Figure  4-1.  The 
scattering  object  lies  at  the  coordinate  origin. 
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E*  ■  x  E(t-z/c) 

Figure  4-1.  The  geometry  of  the  scattering  problem. 

We  will  assume  the  case  of  backscatter  and  r»x.  While  these 
assumptions  are  not  necessary  for  the  general  derivation  of  the  Impulse 
response  they  are  Implied  In  the  use  of  o. 

So  long  as  we  restrict  the  polarization  of  the  Incident  and 
scattered  fields  to  a  single  component,  the  scalar  ratio  E$(jw)/Ef (ju) 
Is  just  a  linear  system  transfer  function.  If  we  define  a  normalized 
transfer  function,  F(jw),  as 

Re  {F(j«)eJ  lt"  /c'}  ■  —  “jTfoj - »  *  *  xcos»(t-z/c),  (4.2) 

then 

o  ■  wc2  |F(Jw)|2  •  U*3) 


27 


The  normalized  transfer  function,  F(J»),  has  a  unique  Inverse 
Fourier  transform  which  Is  the  tine  response  of  the  scatterer.  If  the 
Incident  field  Is  an  Impulse,  6(t-z/c),  the  time  response  Is  the  Impulse 
response,  Fjft).  To  make  sense  physically  we  must  require  that  the 
scattering  problem  be  causal,  that  Is,  there  can  be  no  scattered  field 
before  the  Incident  field  reaches  the  nearest  point  on  the  scatterer. 

In  other  words,  Fj(t)  Is  Identically  equal  to  zero  for  all  t<0,  and  we 
have  the  transform  relation, 

m 

F(Ju)  -  }  F i C t)  e-M  dt  .  (4.4) 

o 

The  time  response  of  a  scatterer  to  an  arbitrary  Input  waveform  may 
be  obtained  by  convolution  with  Fj(t).  Thus  the  Impulse  response 
summarizes  the  response  of  a  scatterer  to  all  Interrogating  waves  at  all 
frequencies.  The  Impulse  response  concept  can  be  applied  to  scatterers 
of  aqy  material  although  we  will  limit  the  discussion  to  perfect 
conductors. 

Some  useful  properties  of  the  Impulse  response  can  be  derived  from 
the  behavior  of  the  scattered  field  In  the  low  frequency  (Rayleigh) 
limit  and  In  the  high  frequency  (physical  optics)  limit.  The  first 
three  moment  conditions  on  Fj(t)  are  obtained  by  enforcing  the  Rayleigh 


law  of  scattering  on  F(J«)  [4]. 

m 

J  Fj(t)  dt  ■  0 

(4.5) 

0 

m 

i  tFj(t)  dt  ■  0 

(4.6) 

o 
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(4  7) 


/  t2Fx(t)  dt  =  2*2 
o 

where  a2  Is  the  coefficient  of  m2  In  an  expansion  of  F(ju>)  near  u>  =  0. 
From  Equation  (4.5)  we  see  that  the  net  area  beneath  the  Impulse 
response  waveform  is  zero. 

In  the  high  frequency  limit  the  current  distribution  over  the 

A 

illuminated  portion  of  the  scatterer  can  be  approximated  by  2n  x  H* . 

The  resulting  physical  optics  approximation  to  the  impulse  response  Is 

F .It)  =  -  1  d2A{z)  .  z  =  ct  ,  (4.8) 

I  w__r_  T- 

where  A(z)  Is  the  cross  sectional  area  Intercepted  by  a  plane  at  z 
parallel  to  the  xy  plane  and  travelling  in  the  positive  z  direction  at  a 
speed  equal  to  one-half  the  velocity  of  the  incident  wave.  A(z)  Is  zero 
for  z<0  and  increases  monotonlcally  to  a  constant  value  at  the  Incident 
shadow  boundary. 

The  physical  optics  approximation  to  Fj ( t)  Is  Insensitive  to  the 
polarization  of  the  Incident  field,  erroneously  predicts  a  discontinuity 
In  Fj ( t )  at  the  shadow  boundary,  and  does  not  predict  any  creeping  wave 
contribution.  But,  In  keeping  with  the  final  value  theorem  of  the 
Laplace  transform,  the  high  frequency  approximation  to  F(jw)  properly 
predicts  the  behavior  of  Fj(t)  near  t  =  0.  Thus  It  provides  a  good 
first  order  approximation  which  can  be  corrected  using  the  moment 
conditions  and  Information  about  creeping  wave  behavior.  A  much  more 
accurate  approximation  may  be  obtained  from  the  Uniform  Theory  of 
Diffraction. 
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Because  of  the  simple  relation  between  the  physical  optics 
approximation  to  Fj ( t)  and  the  geometry  of  the  scattering  object  It  Is 
Instructive  to  derive  the  physical  optics  approximation  for  a  few  simple 
shapes.  A  sphere,  a  right  circular  cylinder,  and  a  disc  represent  three 
types  of  surface  curvature.  The  sphere  has  a  finite  radius  of  curvature 
in  both  principal  planes.  At  broadside  the  cylinder  has  an  infinite 
radius  of  curvature  In  one  principal  plane  and  the  disc  has  Infinite 
radii  of  curvature  In  both.  It  Is  therefore  interesting  to  compare  the 
impulse  response  of  each  near  t  =  0.  The  physical  optics  approximations 
to  the  impulse  response  of  these  three  scatterers  are  derived  In 
Appendix  A  and  illustrated  in  Figure  4-2.  For  convenience  the  time  axis 
is  given  In  units  of  t/r  where  t  equals  the  two  way  free  space  transit 
time  required  to  travel  a  distance  2a,  and  a  is  the  radius  of  the 
sphere,  cylinder,  and  disc. 

For  the  sphere  (Figure  4-2a)  physical  optics  predicts  a  negative 
impulse  of  strength  a/c  followed  by  a  positive  step  of  height  1/2  and 
duration  t.  Correction  to  this  first  order  approximation  will  be 
discussed  shortly. 

Figure  4-2b  gives  the  physical  optics  approximation  to  the  impulse 
response  of  a  2:1  right  circular  cylinder  of  radius  a  at  broadside 
incidence.  At  time  t  =  0  a  negative  impulse  is  followed  by  a  positive 
jump  discontinuity  which  drops  rapidly  to  a  constant  level  of  2/w .  The 

_3 

positive  portion  of  F.(t)  for  0  <  t  <;  r  behaves  as  £  [2tt  -  t2]  ^2 

^  IT 


(  BROADSIDE  ) 


Figure  4-2.  Physical  optics  approximation  to  the  Impulse  response 
of  three  simple  shapes  representing  three  types 
of  surface  curvature  at  the  specular. 
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where  t  =  2a/c.  As  t  approaches  0+  the  singularity  Is  of  the  order 
1  (t)3/2 

T  .  Since  the  positive  excursion  does  not  enclose  a  finite 
/7F 

area  in  zero  time  we  do  not  have  a  true  doublet  t  =  0. 

In  Figure  4-2c  the  physical  optics  approximation  to  Fj ( t )  for  the 

disc  at  broadside  is  shown  to  be  a  doublet  of  strength  (aj2  .  Of 

c 

interest  here  is  the  progression  in  the  overall  shape  of  Fj(t)  due  to 
the  change  in  surface  curvature  at  the  specular  point.  These  three  types 
of  behavior  will  be  exhibited  in  the  measured  impulse  responses. 

Improvements  over  the  physical  optics  approximation  to  Fj(t)  have 
combined  some  rather  ingenious  time  domain  techniques  [3]  with  available 
spectral  solutions  [4],  [7],  The  Geometrical  Theory  of  Diffraction  has 
played  a  key  role.  The  high  frequency  solution  is  able  to  generate 
complex  spectral  data  across  a  wide  enough  band  of  harmonics  without 
excessive  computational  effort.  If  all  significant  scattering  centers 
are  included  in  the  GTD  solution  all  but  the  low  frequency  geometrical 
resonances  can  be  recovered  in  the  impulse  response. 

A  corrected  approximate  impulse  response  for  the  perfectly 
conducting  sphere  is  reproduced  from  [4]  in  Figure  4-3.  The  figure 
compares  the  physical  optics  approximation  with  that  obtained  by  an 
inverse  Fourier  transform  of  475  harmonics  from  an  exact  eigenfunction 
solution.  The  shape  of  the  creeping  wave  may  be  interpreted  from  a  UTD 
analysis.  All  creeping  wave  paths  pass  through  a  point  on  the  shadowed 
side  of  the  sphere  which  is  caustic  for  the  far  field.  The  minimum  path 
gives  the  time  of  the  creeping  wave  peak  at  tcW  =  (2+ir)a/c  =  2.57x.  The 
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Figure  4-3.  Approximate  Impulse  response  of  a  conducting  sphere 
Reproduced  from  [41. 


trailing  portion  Is  time  delayed  due  to  the  phase  factor,  expC-c^tw)], 

found  In  the  UTD  expression  for  the  creeping  wave  [17].  The  precursor 

j*/2 

to  tcw  is  explained  by  the  fact  that  the  UTD  phase  factor,  e  ,  due  to 
crossing  the  caustic  becomes  frequency  dependent  [24]  at  lower 
frequencies  where  UTD  is  no  longer  applicable.  The  combined  effect  Is 
the  characteristic  cusp  shape  exhibited  in  Figure  4-3. 

Before  presenting  measured  impulse  responses  we  must  note  a  shift 
in  the  time  axis  convention  in  the  measured  waveforms.  For  convenience 
in  mounting  and  rotation  the  phase  center  is  taken  to  be  the  center  of 
mass  of  the  scatterer.  This  results  in  a  time  shift  so  that  time  t  *  0 
corresponds  to  the  physical  center  of  the  object.  The  relation  to  the 
time  axis  in  the  theoretical  waveforms  is  obvious  and  should  result  in 
no  confusion. 

The  measured  impulse  response  for  the  sphere  is  presented  in  Figure 
4-4  along  with  an  inverse  transform  of  the  exact  eigenfunction  solution 
covering  the  same  bandwidth  as  the  measurements.  The  negative  specular 
impulse,  the  positive  step,  and  the  creeping  wave  are  all  clearly 
visible.  The  positive  overshoot  at  the  leading  edge  of  the  step  is  a 
Gibb's  phenomenon  due  to  truncation  of  the  spectrum.  The  corresponding 
spectra  are  shown  in  Figures  3-4  and  3-5. 

A  scatterer  which  exhibits  all  three  types  of  surface  curvature  is 
the  right  circular  cylinder  with  a  hemispherical  cap  on  one  end.  The 
geometry  is  shown  in  Figure  2-4.  The  time  axis  of  the  measured  Impulse 
response  waveforms  is  scaled  such  that  x  s  5a/c  is  the  time  required  to 
travel  the  total  length,  5a,  of  the  sphere-capped  cylinder.  The  case  of 
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Figure  4*4.  Measured  Impulse  response  of  a  conducting  sphere 
compared  to  Inverse  transform  of  exact  solution. 
Time  axis  In  units  of  t  ■  2a/c  where  a  Is  sphere 
radius. 

8andw1dth:  383  harmonics  for  0.4  <  ka  <  19. 
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axial  Incidence  on  the  spherical  end  Is  taken  to  be  the  0s  case  while 
the  flat  end  Is  180°. 

The  0°  case  Is  presented  In  Figure  4-5.  Up  to  time  t  *  2a  •  2  T 

c  5 

after  the  specular  return  the  behavior  Is  Identical  to  that  of  the 

sphere.  But  then  the  waveform  drops  rapidly  to  zero  along  the 

cylindrical  sides  until  the  Incident  field  reaches  the  rear  edge.  The 

positive  spike  is  the  direct  diffraction  from  the  rear  edge  and  the 

negative  excursion  Is  the  doubly  diffracted  ray  which  crosses  the  rear 

disc  diameter.  The  peak  of  the  negative  term  occurs  at  a  time  after 

the  positive  spike,  which  Is  the  time  required  for  crossing  one  disc 

jw/4 

diameter.  The  doubly  diffracted  ray  picks  up  a  phase  term  e  from 

jw/2 

each  of  the  two  diffractions  and  another  e  in  crossing  the  caustic 

on  the  cylinder  axis.  This  gives  it  a  polarity  opposite  to  that  of  the 

single  diffraction.  With  respect  to  the  time  resolution  shown  In  Figure 

3-2,  note  that  the  singly  diffracted  component  Is  of  minimum  width, 

indicating  It  has  been  diffracted  as  If  from  a  single  stationary  point, 

while  the  doubly  diffracted  component  exhibits  the  time  dispersion 

connected  with  crossing  the  caustic  on  the  cylinder  axis. 

The  measured  Impulse  response  of  the  sphere-capped  cylinder  may  be 

compared  to  a  theoretical  approximation  obtained  by  a  space-time 

♦ 

Integral  equation  approach  [8],  In  this  technique  the  H-fleld  Integral 
equation  Is  solved  In  the  time  domain  for  the  surface  currents  on  the 
scattering  object,  and  the  currents  then  used  to  compute  the  scattered 
field.  An  Impulse  response  augmentation  technique  [8]  can  then  be  used 
to  Improve  the  quality  of  the  singular  (high  frequency)  portions  of  the 
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Figure  4-5.  Measured  Impulse  response  of  a  conducting  2:1  right 
circular  cylinder  with  hemispherical  cap  on  one  end 
0°  case:  axial  Incidence  on  the  sphere-capped  end. 
Time  axis  In  units  of  t  *  ®a/c  where  L  *  5a  Is  the 
total  length. 

Bandwidth:  960  harmonics  for  1  <  kL  <  96. 


scattered  waveform.  The  resulting  Impulse  response  for  axial  Incidence 
on  the  sphere-capped  end  Is  shown  In  Figure  4-6a.  The  comparison  with 
the  measured  results  In  Figure  4-5  Is  quite  favorable  except  for  the 
relative  magnitude  of  the  two  diffracted  components  from  the  rear  edge. 
Again  note  the  difference  In  the  widths  of  these  two  components. 

Turning  now  to  axial  Incidence  on  the  sphere-capped  cylinder  at 
180°,  the  specular  surface  is  a  flat  disc.  The  measured  Impulse 
response  at  180°  is  presented  In  Figure  4-7.  As  expected  from  the 
physical  optics  approximation  the  specular  return  from  the  flat  end  of 
the  cylinder  Is  a  doublet.  But  following  the  specular  and  ringing  for  a 
duration  of  at  least  2t  appears  a  wave  with  a  resonant  time  between 
positive  and  negative  peaks  of  2a/Cf  corresponding  to  one  disc  diameter. 
Even  the  creeping  wave,  which  we  at  expect  at  a  time  tcW  =  (10  +  w)a/c  = 
2.6t  after  the  specular,  appears  distorted  by  this  ringing.  We  would 
expect  the  creeping  wave  to  be  comparable  In  shape  and  magnitude  to  that 
from  a  sphere  of  the  same  radius  (shown  In  the  inset),  since  the 
cylinder  walls  should  have  little  effect  on  the  field. 

An  equivalent  currents  solution  [9]  for  the  backscatter  from  a 
perfectly  conducting  circular  disc  of  radius  a  at  normal  Incidence  gives 
a  scattered  field  transfer  function, 

1-1  -jf2a)s 

F(  s)=  [if  S  +  i_  (l)7  s  7  e  c  ,  s  =  ju,  .  (4.9) 

c  /it  c 

Using  the  Laplace  transform  pair, 

1  - *  JL  •  (4.10) 

/irt  /S 
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Figure  4-6.  Theoretical  Impulse  response  approximation  obtained 
via  space-time  Integral  equation  technique  [81  for  a 
conducting  2:1  right  circular  cyllhder  with 
hemispherical  cap  on  one  end,  at  axial  Incidence. 
Bandwidth:  0  <  ka  <  30. 
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Figure  4-7.  Measured  Impulse  response  of  a  conducting  2:1  right 
circular  cylinder  with  hemispherical  cap  on  one  end 
180°  case:  axial  Incidence  on  the  flat  end. 

Time  axis  In  units  of  t  ■  5a /c  where  L  -  5a  Is  the 
total  length. 

Bandwidth:  960  harmonics  for  1  <  kL  <  96. 


and  the  time  shift  theorem  we  find  the  inverse  transform  of  F(s)  given 

by 

1 

Ft  (t )  *  fi.)2  «'(t)  +1  fa)7  1  »  f4-11) 

1  c  »  c  /t  -  5i 7" 

The  first  term  is  just  the  physical  optics  prediction.  The  second  term 
suggests  a  resonant  peak  at  t  =  2a/c. 

In  [10]  the  exact  Wiener-Hopf  solution  for  scattering  by  an  open 
circular  waveguide  in  the  resonance  region  is  transformed  into  a  time 
domain  expression  for  the  scattered  electric  field  on  axis  given  by 
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(4.12) 


Since  the  cylinder  Is  open  ended  we  would  not  expect  the  first  term  to 
be  a  doublet  as  It  Is  for  the  flat  end  of  the  cylinder.  But  the  second 
term,  which  has  resonant  peaks  at  multiples  of  2a/c,  appears  to  be 
associated  with  multiple  diffraction  across  the  diameter.  In  Figure  4-8 
the  normalized  time  domain  waveform  is  reproduced  from  [10].  The 
authors  have  noted  that  the  ringing  becomes  negligible  after  the  third 
peak.  The  same  is  true  for  the  measured  Impulse  response  of  Figure  4-7. 
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Finally,  returning  to  the  space-time  Integral  equation  approach, 
from  [8]  the  Impulse  response  for  the  sphere-capped  cylinder  with  axial 
Incidence  on  the  flat  end  Is  shown  In  Figure  4-6b.  Here  the  resonance 
Is  not  as  pronounced  and  does  not  continue  through  the  creeping  wave. 
Some  of  the  ringing  after  the  Initial  doublet  can  be  seen,  and  the 
rounded  positive  swing  preceedlng  the  creeping  wave  Is  very  similar  to 
the  corresponding  portion  of  the  measured  Impulse  response. 

So  far  we  have  compared  measured  time  waveforms  with  several 
theoretical  impulse  responses.  All  are  approximate  in  that  they 
represent  discrete  Fourier  transforms  of  bandllmlted  spectra.  The 
measured  spectra.  In  addition,  are  corrupted  by  noise  which  has  been 
extracted  to  some  extent  by  the  methods  described  in  Chapter  III.  The 
Gibb's  phenomenon  associated  with  truncation  of  the  spectrum  and  other 
discontinuities  has  been  treated  to  varying  degrees  in  all  the  time 
waveforms.  When  viewing  the  measured  waveforms  It  should  be  noted  that 
the  residual  high  frequency  ripple  may  be  further  suppressed  at  the 
expense  of  time  resolution. 

For  the  sphere  and  for  axial  Incidence  on  a  body  of  revolution  the 
backscattered  radar  cross  section  Is  Independent  of  the  polarization  of 
the  Incident  field.  When  the  radar  moves  off  axis  the  two  polarization 
cases  must  be  treated  separately.  The  case  of  the  Incident  electric 
field  perpendicular  to  the  plane  formed  by  the  direction  of  Incidence 
and  the  axis  of  the  scatterer  Is  designated  perpendicular  (TE) 
polarization.  The  case  of  the  Incident  electric  field  lying  In  that 
plnae  Is  designated  parallel  (TM)  polarization. 


The  measured  impulse  response  for  the  sphere-capped  cylinder  at 

broadside  (90°)  is  shown  in  Figure  4-9  for  perpendicular  (TE) 

polarization  and  Figure  4-10  for  parallel  (TM)  polarization.  The 

analysis  splits  most  easily  into  two  parts:  behavior  near  the  specular 

and  behavior  for  late  time.  Near  the  specular  the  waveform  should 

behave  much  like  the  physical  optics  approximation  for  the  right 

circular  cylinder  without  the  spherical  cap  given  In  Figure  4-2b,  where 

we  noted  a  negative  impulse  and  a  positive  singularity  of  order  (£)3/2 

c 

tf3/2.  A  consideration  of  the  surface  currents  on  the  shadowed  side  of 
the  cylinder  [11]  would  lead  us  to  expect  physical  optics  to  approximate 
the  TM  case  more  closely  than  the  TE  case.  In  the  TM  waveform  of  Figure 
4-10  we  see  a  negative  Impulse  at  the  specular  and  a  positive 
singularity  less  than  half  as  high  which  drops  rapidly  to  a  positive 
constant  value  until  a  time  corresponding  to  a  distance  of  one  cylinder 
diameter  after  the  specular.  This  Is  the  same  overall  shape  seen  In  the 
physical  optics  approximation. 

A  first  order  Uniform  Theory  of  Diffraction  solution  for  the 
sphere-capped  cylinder  [12]  has  been  inversed  Fourier  transformed  to 
produce  the  Impulse  response  approximation  In  Figure  4-11.  The  geometry 
for  the  first  order  UTD  solution  Is  shown  In  Figure  A-4  for  the  various 
angles  of  incidence.  The  first  order  solution  Includes  the  contribution 
from  the  edge  at  the  flat  end  of  the  cylinder  and  the  discontinuity  In 
surface  curvature  at  the  junction  of  the  spherical  cap.  At  aspects  off 
broadside  a  third  contribution  comes  from  either  the  reflected  specular 
on  the  sphere  or  the  rear  edge  of  the  flat  end  rim  depending  upon  which 
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Figure  4-9.  Measured  Impulse  response  of  a  conducting  2:1  right 
circular  cylinder  with  hemispherical  cap  on  one  end. 
90°  case:  broadside. 

Perpendicular  polarization  (TE). 

Time  axis  in  units  of  t  >  5a/c  where  l»5a  is  the  total 
length. 

Bandwidth:  960  harmonics  for  1  <  kL  <  96. 


Measured  Impulse  response  of  a  conducting  2:1  right 
circular  cylinder  with  hemispherical  cap  on  one  end. 
90°  case:  broadside. 

Parallel  polarization  (TM). 

Time  axis  In  units  of  t  *  *a/c  where  L*5a  Is  the  total 
length. 

Bandwidth:  960  harmonics  for  1  <  kL  <  96. 


First  order  UTD  Impulse  response  of  a  conducting  2:1 
right  circular  cylinder  with  hemispherical  cap  on  one  end 
90°  case:  broadside. 

Parallel  polarization  (TM), 

Time  axis  In  units  of  t  «  5a/c  where  l*5a  is  the  total 
length. 

Bandwidth:  960  harmonics  for  2  <  kL  <  96. 


lies  in  the  lit  region.  A  thorough  discussion  of  the  UTD  solution  Is 
found  in  [12].  In  the  symmetrical  case  of  the  cylinder  with  spherical 
caps  on  both  ends  the  transfer  function  for  the  scattered  field  at 
broadside  simplifies  to 

F(jw)  =  -  f8l^/2  /juj  .  (4.13) 

C  /i  C 

We  can  see  that  the  first  term  gives  the  Inverse  transform  .  (aj  6(t)» 

c 

which  is  the  specular  from  a  sphere.  The  term  Involving  /3u>,  associated 
with  the  cylindrical  curvature,  does  not  have  a  well  defined  inverse 
Laplace  transform  In  analytic  form  because  of  a  singularity  of  order 
t"^/2.  But  this  is  exactly  the  behavior  we  have  noted  previously  for 
times  immediately  following  the  specular  return. 

Because  of  the  high  quality  of  the  UTD  spectrum  the  amplitude  of 
the  UTD  impulse  response  in  Figure  4-11  has  been  reduced  to  3/4  scale. 
Allowing  for  this  scale  factor  the  measured  TM  impulse  response  compares 
extremely  well  with  the  UTD  waveform. 

Continuing  the  early  time  analysis,  let  us  consider  the  measured 
TE  impulse  response  in  Figure  4-9.  The  negative  impulse  for  this 
polarization  Is  nearly  identical  to  that  of  the  TM  case.  But  the 
positive  singularity  has  significantly  greater  amplitude  and  drops  to  a 
positive  step  level  much  higher  than  the  TM  waveform.  To  understand 
this  behavior,  refer  to  [7]  In  which  the  impulse  response  of  a  right 
circular  cylinder  Is  derived  from  a  normalization  of  the  Keller  GTD 
solution  for  the  Infinite  cylinder  [13].  The  results,  shown  in  Figure 
4-13,  give  a  sign  difference  between  the  two  polarizations  in  the  jump 
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discontinuity  after  the  specular.  The  difference  is  based  on  the  sign 
difference  in  the  second  term  of  the  Luneberg-Kl 1 ne  expansion  for  the 
reflected  field.  This  term,  involving  powers  of  1/jw,  Is  a  higher  order 
term  which  Is  not  Included  in  the  first  order  UTD  solution.  Its  inverse 
transform  is  a  unit  step,  u(tl.  For  the  infinite  cylinder  the  reflected 
field  is  the  only  contribution  to  the  scattered  field  except  for 
possible  creeping  waves  which  would  appear  later  in  time.  But  when  the 
cylinder  is  finite  the  end  contributions,  which  we  have  seen  to  be  of 
order  /Ju  dominate  the  l/jw  term.  It  Is  just  this  l/ju>  term  which 
accounts  for  the  difference  between  the  two  measured  impulse  responses 
at  the  jump  discontinuity  Immediately  after  the  specular  return. 

The  conclusions  to  be  drawn  at  this  point  are  twofold.  First,  the 
measurements  have  recovered  a  higher  order  effect.  Secondly,  while  the 
first  order  UTD  solution  converges  at  broadside  to  the  same  results  for 
the  TE  and  TM  cases,  a  UTD  solution  whch  Includes  the  higher  order  term 
is  quite  capable  of  predicting  the  slight  difference  between  the  two 
polarizations. 

The  late  time  analysis  will  begin  with  the  measured  7W  impulse 
response  in  Figure  4-10.  The  waveform  reaches  a  relative  null  at  a  time 
corresponding  to  approximately  one  half  the  total  object  length.  At 
this  point  a  low  frequency  whole-body  resonance  begins  which  dominates 
the  remainder  of  the  waveform.  The  first  peak  of  this  resonance  Is  at 
approximately  4/5t  after  the  specular  and  the  resonant  period  is 
slightly  less  than  0.7t.  a  complete  solution  for  the  scattered  field 
from  a  cylindrical  structure  with  length  to  diameter  ratio  on  the  order 
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of  2:1  is  not  yet  available  in  the  resonance  region.  But  as  a  rough 
approximation  the  sphere-capped  cylinder  may  be  modeled  as  a  short  fat 
dipole.  Calculations  of  the  radar  cross  section  of  cylindrical  antennas 
are  given  in  [14]  as  a  function  of  frequency  and  length  to  radius  ratio 
(See  also  [15]).  While  these  results  are  given  only  for  length  to 
radius  ratios  as  low  as  L/a  =  45  (compared  to  L/a  *  5  for  the 
sphere-capped  cylinder)  the  general  trend  indicates  a  second  resonant 
peak  at  kL  =  3n  which  becomes  much  stronger  for  decreasing  L/a  ratios. 

As  can  be  seen  in  Figure  4-10,  the  measured  waveform  resonates  with  a 
period  T2  =  0.7r,  or  kL  =  2.9n. 

The  entire  TM  waveform  appears  to  be  superimposed  upon  a  sinusoid 

of  period,  T^,  equal  to  approximately  3t,  corresponding  to  kL  * 

While  the  first  resonant  peak  for  a  thin  dipole  occurs  at  kL=n  (or  L  = 
x/2)  it  is  noted  In  [14]  that  this  peak  shifts  noticeably  toward  lower 
values  of  kL  as  the  length  to  radius  ratio  decreases.  This  Is 
consistent  with  the  period,  Tj,  noted  above.  Thus  the  measured  TM 
waveform  shows  evidence  of  two  natural  resonances  of  the  sphere-capped 
cylinder. 

By  means  of  comparison  the  measured  TM  Impulse  response  of  a  2:1 
right  circular  cylinder  without  the  spherical  cap  is  shown  in  Figure 
4-12.  All  the  characteristics  of  the  TM  sphere-capped  cylinder  waveform 
are  again  evident  in  Figure  4-12,  except  that  the  natural  resonance 
appears  to  be  not  quite  as  strong.  We  would  expect  this  to  be  the  case 
since  the  spherical  cap  decreases  the  Influence  of  the  radial  currents 
across  the  ends  of  the  cylinder.  For  this  same  reason  the  dipole  model 
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Figure  4-12.  Measured  Impulse  response  of  a  conducting  2:1  right 
circular  cylinder. 

90°  case:  broadside. 

Parallel  polarization  (TM). 

Time  axis  In  units  of  t  -  **/c  where  L»4a  Is  the  total 
length. 

Bandwidth:  766  harmonics  for  0.8  <  kL  <  76. 


more  closely  approximates  the  sphere-capped  cylinder  than  the  flat  end 
cylinder  [16]. 

The  late  time  behavior  of  the  measured  TE  impulse  response  of  the 
sphere-capped  cylinder  in  Figure  4-9  has  a  positive  step  and  a  creeping 
wave.  The  creeping  wave  around  the  cylinder  appears  in  Figure  4-13, 
which  is  reproduced  from  [7].  In  Figure  4-13  the  specular  term  has  been 
removed  and  the  polarity  of  the  perpendicular  (TE)  waveform  reversed. 
With  the  specular  term  removed  It  can  be  seen  that  the  second  term  of 
the  Luneberg-Kl ine  expansion  is  responsible  for  the  step  behavior.  In 
both  Figures  4-9  and  4-13  the  transition  from  the  positive  step  into  the 
creeping  wave  appears  more  abrupt  for  the  cylindrical  surface  than  for 
the  sphere,  due  to  the  fact  that  the  creeping  wave  paths  do  not  cross  a 
caustic. 

Superimposed  upon  the  step  portion  of  the  measured  waveform  of 
Figure  4-9  is  a  pronounced  ringing.  In  light  of  the  amount  of  precursor 
In  this  particular  waveform  it  appears  that  the  ringing  is  measurement 
error.  In  this  connection  the  measured  TE  Impulse  response  for  the 
right  circular  cylinder  without  the  spherical  end  cap  is  shown  In  Figure 
4-14.  The  aforementioned  ringing  Is  noticeably  absent.  All  other 
features  describes  for  the  sphere-capped  cylinder  are  in  evidence.  Note 
in  particular  the  rapid  transition  Into  the  creeping  wave. 

When  the  sphere-capped  cylinder  is  rotated  to  a  45°  look  angle  the 
scattering  phenomena  become  slightly  more  complicated.  The  measured  TM 
Impulse  response  at  45°  Is  shown  In  Figure  4-15  along  with  the 
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Figure  4-13.  GTD  Impulse  response  of  a  conducting  right  circular 
cylinder  (reproduced  from  [7]). 

Specular  term  removed. 

90°  case:  broadside. 

Polarity  of  perpendicular  polarization  waveform  reversed. 
Time  axis  In  units  of  t  ■  2a/c  where  a  ■  cylinder 
radius. 

Bandwidth:  79  harmonics  for  0.157  <  ka  <  12.4. 
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Figure  4-14.  Measured  Impulse  response  of  a  conducting  2:1  right 
circular  cylinder. 

90°  case:  broadside. 

Perpendicular  polarization  (TE). 

Time  axis  in  units  of  t  «  *a/c  where  L*4a  is  the  total 
length. 

Bandwidth:  766  harmonics  for  0.8  <  kl  <  76. 


corresponding  first  order  UTD  Impulse  response.  The  TE  case  Is  shown  In 
Figure  4-16.  Note  that  the  measured  TM  impulse  response  covers  a 
smaller  bandwidth  (  1.0  <  kL  <  64.  )  than  the  TE  case.  In  order  to 
obtain  a  fair  comparison  the  UTD  solution  is  restricted  to  the  same 
bandwidth.  Again  the  analysis  is  split  Into  times  near  the  specular  and 
late  times. 

The  specular  return  Is  identical  for  both  polarizations  and  for  UTD 
versus  measured  waveforms,  indicating  the  specular  term  is  polarization 
independent.  Immediately  after  the  specular  both  measured  waveforms 
begin  a  positive  step  just  as  In  the  case  of  the  sphere  or  sphere-capped 
cylinder  at  0°.  This  step  behavior  is  not  seen  In  the  first  order  UTD 
waveform  indicating  the  absence  of  the  1/k  term  from  the  Luneberg-Kl ine 
expansion,  as  discussed  In  connection  with  the  broadside  cylinder. 

When  the  Incident  field  reaches  the  junction  between  the  sphere  and 
cylinder  we  see  a  noticeable  difference  In  the  response  of  the  two 
polarizations  to  the  discontinuity  In  surface  curvature  at  the  junction. 
As  a  heuristic  explanation  consider  that  in  the  TE  case  the  radius  of 
curvature,  a  ,  in  the  direction  parallel  to  the  incident  electric  field 
vector  does  not  change  across  the  junction  between  sphere  and  cylinder. 
The  TE  waveform  responds  relatively  slowly  and  continuously  to  the 
discontinuity  in  surface  curvature.  In  the  TM  case  the  radius  of 
curvature  In  the  direction  parallel  to  the  Incident  electric  field 
vector  jumps  from  a  to  Infinity  across  the  junction.  In  the  measured  TM 
waveform  the  response  drops  abrubtly  to  zero  at  the  Junction  Indicating 
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Figure  4-15.  Measured  versus  first  order  UTD  Impulse  response  of  a 

conducting  2:1  right  circular  cylinder  with  hemispherical 
cap  on  one  end. 

45°  case. 

Parallel  polarization  (TM), 

Time  axis  In  units  of  t  ■  ®a/c  where  L«5a  Is  the  total 
length. 

Bandwidth:  640  harmonics  for  1  <  kL  <  64. 
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Figure  4-16. 


Measured  versus  first  order  UTD  Impulse  response  of  a 
conducting  2:1  right  circular  cylinder  with  hemispherical 
cap  on  one  end. 

45°  case. 

Perpendicular  polarization  (TE). 

Time  axis  In  units  of  t  ■  5a/c  where  L*5a  Is  the  total 
length. 

Bandwidth:  960  harmonics  for  1  <  kL  <  96. 


the  Impulse  response  fs  polarization  sensitive  to  the  discontinuity  In 
surface  curvature. 

The  late  time  behavior  shows  significant  differences  between  the 
two  polarizations  associated  with  diffraction  from  the  rear  edge  of  the 
cylinder.  We  will  begin  the  late  time  analysis  with  the  TE  case.  The 
UTD  waveform  shows  a  weak  positive  return  from  the  point  of  the  rear 
cylinder  rim  nearest  the  radar.  This  small  spike  appears  as  the  high 
ooint  in  the  center  of  a  wide,  shallow  positive  swing  in  the  measured 
waveform.  The  mechanism  responsible  appears  to  be  a  low  frequency 
effect  outside  the  regime  of  UTD  and  localized  scattering. 

The  distinctive  late  time  feature  in  the  measured  TE  waveform  is  a 
wide  negative  excursion  centered  at  a  time  2t  after  the  specular.  It 
has  the  time  dispersion  and  cusp  shape  characteristic  of  a  creeping 
wave.  A  plausible  explanation  for  this  scattering  mechanism  may  be 
proposed  using  a  UTD  argument  involving  a  spiral  creeping  wave  path  and 
diffraction  from  the  rear  edge.  Figure  A-17  shows  the  geometry  for  the 
possible  ray  paths.  The  incident  shadow  boundary  lies  along  the  top  of 
the  cylinder  parallel  to  the  axis.  A  continuum  of  rays  may  attach  at 
this  shadow  boundary  and  creep  along  a  spiral  path  to  the  rear  cylinder 
edge.  Each  ray  diffracts  from  the  edge  and  creeps  back  along  a 
symmetric  spiral  path  on  the  other  side  of  the  cylinder  to  launch  in  a 
backscatter  direction  at  the  opposite  shadow  boundary.  If  the  cylinder 
is  "unwrapped"  from  shadow  boundary  to  shadow  boundary  it  can  be  shown 
by  simple  geometry  that  the  same  total  distance  is  travelled  by  each 
of  these  rays.  The  path  length  referenced  to  the  specular  point  is  i  * 
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2a [  1+2/7  +  *!'in\  -  9.88a.  Thus  we  would  expect  the  peak  of  the 
creeping  wave  to  occur  at  a  time  ^  =  1.98t  after  the  specular.  Due  to 

1  /3 

the  phase  factor  e"^*ctcw"^w  ^  discussed  earlier  in  connection  with 
creeping  waves  [17],  we  expect  this  feature  to  he  dispersed  in  time. 

The  difference  between  the  expected  time  of  arrival  of  this  proposed 
mechanism  and  the  time  read  off  the  measured  waveform  is  less  than  1 
percent.  Such  a  comparison  is  very  convincing  evidence  that  the 
proposed  mechanism  accounts  for  the  measured  phenomenon. 

The  foregoing  discussion  may  serve  to  demonstrate  the  utility  of 
the  impulse  response  waveform  in  understanding  specific  scattering 
mechanisms  and  improving  the  accuracy  of  theoretical  formulations.  The 
expected  arrival  time  of  the  response  from  a  specific  scattering  center 
is  obtained  from  exactly  the  same  calculation  as  the  phase  factor 
associated  with  each  diffraction  coefficient  in  a  UTD  analysis. 

Therefore  a  logical  and  convenient  place  to  verify  the  UTD  results  is  on 
the  time  axis.  Since  UTD  is  a  superposition  of  distinct,  local 
scattering  mechanisms  the  impulse  response  displays  their  relation  to 
the  geometry  of  the  scatterer  quite  naturally. 

The  late  time  behavior  of  the  TM  waveform  presents  an  equally 
intriguing  problem.  At  the  leading  edge  of  the  rear  cylinder  rim  both 
the  UTD  and  measured  waveforms  show  a  strong  positive  spike.  This  point 
is  the  stationary  point  for  first  order  diffraction  from  the  entire 
rear  rim  of  the  cylinder.  The  measured  waveform  shows  this  single 
diffraction  mechanism  at  the  leading  edge  of  a  wide  positive  excursion. 
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The  position,  height,  and  rise  time  of  this  feature  is  nearly  identical 
to  that  of  the  UTD  waveform.  But  after  peaking  at  the  leading  edge  of 
the  rear  rim  the  measured  TM  waveform  tapers  slowly  and  continuously  to 
zero  after  a  time  of  approximately  0.7t  and  remains  negative  for  another 
t  seconds. 

The  measured  waveform  suggests  that  a  ring  current  is  induced  on 
the  cylinder  rim.  Such  a  mechanism  could  be  treated  with  an  equivalent 
currents  approach  [9].  As  a  step  toward  understanding  the  phenomenon 
the  UTD  impulse  response  has  been  subtracted  from  the  measured  TM 
impulse  response.  The  difference  waveform  is  shown  in  Figure  4-18  along 
with  the  UTD,  measured,  and  difference  spectra.  The  spectrum  of  the 
difference  shows  significant  energy  in  both  the  low  and  high  frequency 
regions.  In  addition,  the  difference  spectrum  shows  a  fairly  even 
periodicity  at  intervals  of  kL  =  5  corresponding  to  a  fundamental 
frequency  of  u0  =  5/t.  It  can  be  seen  that  this  resonance  does  not 
correspond  to  a  period  evident  in  either  the  measured  or  UTD  spectrum. 

If  we  related  the  interference  to  a  path  length,  i,  corresponding  to 
some  physical  dimension  of  the  scatter  we  would  get  nulls  at  intervals 
of  e"Jki  =  e-Jnn.  Using  the  value  of  u>0  we  obtain  i  =  nna  or  half  the 
circumference  of  the  cylinder  rim.  It  would  seem  logical  to  expect  the 
resulting  time  response  to  change  polarity  at  n  =  1  corresponding  to  a 

time  =  0.63t  after  the  leading  edge  of  the  rear  rim.  This  Is  the 
approximate  time  of  the  zero  crossing  in  the  measured  waveform. 
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Figure  4-18.  (a)  Time  waveform  of  difference  between  measured  and 

first  order  UTD  Impulse  responses  for  sphere- 
capped  cylinder  at  45°;  TM. 

(b)  UTD,  measured,  and  difference  spectra  of  sphere- 
capped  cylinder  at  45°;  TM. 
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The  above  discussion  is  somewhat  heuristic  and  is  not  meant  to 
provide  the  definitive  explanation  for  the  observed  behavior  but  is 


Intended  to  illustrate  the  possible  uses  of  the  impulse  response 
waveform  in  scattering  analysis. 

Finally  the  sphere-capped  cylinder  is  rotated  to  135°.  At  this 
look  angle  we  do  not  have  a  reflected  field  specular  as  in  the  previous 
cases.  The  TM  Impulse  response  Is  presented  in  Figure  4-19  and  the  TE 
in  Figure  4-20.  We  will  begin  the  early  time  analysis  with  the  TM 
waveform. 

In  Figure  4-19  we  see  an  excellent  correlation  between  the  measured 
and  UTD  waveforms  at  the  front  cylinder  disc.  The  response  is  very 
similar  to  the  tangent  curve  obtained  from  the  physical  optics 
approximation  for  the  disc  away  from  normal  Incidence  (Figure  B-3).  For 
this  polarization  the  response  is  symmetric  from  both  leading  and 
trailing  edges  of  the  disc.  We  have  already  seen  that  as  the  angle  of 
incidence  rotates  to  normal  the  response  converges  to  a  doublet. 

For  TE  polarization  In  Figure  4-20  we  do  not  see  the  same  symmetry. 
The  response  from  the  leading  edge  is  strong  and  that  from  the  trailing 
edge  much  weaker.  Again  the  UTD  and  measured  waveforms  compare  very 
well.  Note  that  just  as  In  the  broadside  cylinder  case  the  physical 
optics  approximation  provides  a  better  prediction  of  the  TM  behavior. 

We  may  safely  say  that  the  early  time  response  Is  sensitive  to 
polarization. 

While  the  early  time  response  characterizes  the  sphere-capped 
cylinder  for  this  look  angle  the  late  time  response  shows  some  more 
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Measured  versus  first  order  UTD  impulse  response  of  a 
conducting  2:1  right  circular  cylinder  with  hemispherical 
cap  on  one  end. 

135°  case. 

Parallel  polarization  (TM), 

Time  axis  In  units  of  t  ■  ®a/c  where  L»5a  is  the  total 
length. 

Bandwidth:  640  harmonics  for  1  <  kL  <64. 
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Measured  versus  first  order  UTD  Impulse  response  of  a 
conducting  2:1  right  circular  cylinder  with  hemispherical 
cap  on  one  end. 

135°  case. 

Perpendicular  polarization  (TE). 

Time  axis  In  units  of  t  *  5a/c  where  L*5a  Is  the  total 
length. 

Bandwidth:  640  harmonics  for  1  <  kL  <  64. 


subtle  features.  We  do  not  see  any  of  the  resonance  associated  with  the 
disc  diameter  as  we  saw  In  the  180°  case.  But  the  TM  waveform  shows  a 
low  frequency  ringing  with  a  period  of  approximately  0.9t.  We  would 
expect  a  weaker  version  of  the  dipole  resonance  noted  In  the  broadside 
case,  where  the  period  of  0.7t  agreed  with  theoretical  results.  The 
nature  of  such  ringing  for  aspects  off  broadside  Is  not  obvious. 

In  the  TM  waveform  we  see  a  creeping  wave  which  attaches  on  the 
spherical  cap,  creeps  around  the  sphere,  skirts  the  cylinder  wall,  and 
diffracts  from  the  cylinder  edge  furthest  from  the  radar.  From  the 
geometry  we  can  calculate  the  arrival  time  of  this  creeping  wave  as 
tcw  =  1.4t.  This  time  corresponds  to  the  lowest  point  on  the  late  time 
TM  waveform.  The  cusp  shape  Is  apparent  even  though  it  appears 
superimposed  upon  the  low  frequency  sinusoid  noted  earlier.  Although 
the  magnitude  of  this  creeping  wave  is  decreased  by  the  edge  diffraction 
It  Is  still  seen  to  be  a  measurable  contribution. 

The  distinctive  feature  in  the  late  time  TE  waveform  is  the 
creeping  wave  which  attaches  at  the  sphere-cylinder  junction,  creeps 
around  the  sphere,  and  launches  from  the  point  on  the  junction 
diametrically  opposite.  It  is  seen  at  tcW  =  l.lx  In  Figure  4-20.  Its 
shape  and  magnitude  indicates  that  It  Is  unaffected  by  the  presence  of 
the  cylinder.  For  both  polarizations  It  Is  a  simple  matter  to  Include 
the  creeping  wave  term  In  the  UTD  solution. 

The  preceeding  discussion  has  covered  five  major  look  angles  and 
two  principal  polarizations.  The  sphere-capped  cylinder  Is  seen  to 
represent  several  types  of  specular  behavior,  multiple  Interactions,  and 


low  frequency  resonance.  Specific  scattering  mechanisms  may  be  verified 
by  rotating  the  object  slightly  and  examining  the  resultant  time  shift, 
or  absence  of  it,  in  the  Impulse  response.  Multiple  Interactions 
between  the  sphere  and  cylinder  fall  out  when  the  sphere-capped  cylinder 
waveform  Is  compared  to  that  of  the  cylinder  without  the  spherical  cap. 
And  by  subtracting  known  UTD  terms  from  the  measured  waveforms  we  find 
out  which  additional  scattering  mechanisms  are  significant.  These  and 
other  mechanisms  may  be  examined  further  by  considering  the  cross 
polarized  Impulse  response.  All  of  these  methods  are  made  possible  by 
the  unique  relation  between  the  Impulse  response  waveform  and  the 
scattering  geometry. 


CHAPTER  V 


CROSS -POLAR  I ZED  MEASUREMENTS 

In  order  to  measure  cross -polarized  backscatter,  It  was  necessary 
to  select  a  suitable  calibration  target  to  replace  the  sphere,  which  has 
no  cross-polarized  return.  A  calibration  target  must  satisfy  the 
following  criteria:  1)  it  must  be  broadband,  l.e.,  not  highly  resonant; 
2)  it  must  have  a  known,  readily  available  exact  solution;  and  3)  it 
must  be  relatively  insensitive  to  positioning.  The  geometry  for  the 
cross-polarized  measurements  Is  sketched  In  Figure  5-1.  A  narrow,  thin 
strip  at  edge-on  incidence  was  chosen  as  a  calibration  target.  At 
edge-on,  the  thin  strip  was  known  to  satlfy  the  third  criterion  as  shown 
In  Figure  5-2  [21].  It  was  predicted  and  proven  that  a  narrow  strip 
would  be  broadband.  And,  finally,  a  moment  method  solution  was 
available. 

At  a  tilt  of  45°  off  vertical,  the  edge-on  strip  gives  a  cross- 
polarized  return  comparable  In  magnitude  to  the  co-polarlzed  return  of 
the  calibration  spheres.  The  strip  has  a  cosine  sensitivity  to  error  in 
the  amount  of  tilt.  Co-polarized  measurements  of  the  edge-on  strip 
calibrated  against  a  sphere,  shown  In  Figure  5-3,  verify  that  the  narrow 
(<  A/4)  strip  Is  indeed  broadband  In  both  amplitude  and  phase. 
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Figure  5-1.  Calibration  target  for  cross  polarized  measurements 


69 


,3 


s$tn. 


1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

so.  60.  90. 

PH  =  180  —  <f>' 


Figure  5-2.  H-plane  echo  width  pattern  of  X/4  wide  strip. 
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Figure  5-3.  Co-pol arl zed  measurements  of  a  narrow  strip 
calibrated  against  a  sphere. 
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The  targets  described  in  Chapter  II  were  tilted  out  of  the 
principal  planes  in  order  to  ensure  a  cross-polarized  return.  In 
addition  to  these  targets,  several  sizes  of  discs,  flat  plates  and 
strips  were  included.  To  test  the  viability  of  obtaining  cross- 
polarized  results  using  the  strip,  the  measured  spectra  were  compared  to 
the  few  existing  cross-polarized  solutions.  In  Figure  5-4  an  exact 
solution  by  D.  Hodge  [22]  is  compared  to  measured  cross-polarized 
results  for  a  disc  at  edge-on  incidence  and  at  a  45°  tilt.  The  results 
compare  very  well  for  those  portion  of  the  spectrum  above  the  noise 
floor  of  the  system.  While  the  cross-polarized  target  cross  section 
tends  to  be  lower  than  in  the  co-polarized  case,  the  noise  floor  is 
generally  lower  for  the  cross-polarized  measurements  due  to  the 
decoupling  of  the  transmit  and  receive  antennas.  The  cross-polarized 
measurements  show  a  reasonable  degree  of  accuracy  down  to  cross  sections 
of  1  sq.  cm. 

While  the  time  domain  analysis  Is  not  yet  complete,  these  results 
demonstrate  that  the  measurement  system  described  in  Chapter  II  has  been 
successfully  modified  to  measure  cross-polarizced  radar  cross  sections. 
The  edge-on  strip  has  been  shown  to  be  an  accurate  calibration  target, 
and  has  been  used  to  obtain  the  cross-polarized  Impulse  responses  given 
In  [23],  The  cross-polarization  data  complete  the  polarization 
scattering  matrix,  and  in  conjunction  with  the  co-polarized  Impulse 
responses,  completely  characterized  the  scattering  behavior  of  the 
targets  in  question. 
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17.7  CM  DIflM.  DISC  EDGE  ON  45DEG  TILT  2-4GHZ  CP 


gure  5-4.  Comparison  of  measured  and  theoretical  (Hodge)  cross 
polarized  response  from  a  thin  disk  tilted  at  45°. 
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CHAPTER  VI 


CONCLUSIONS 

The  impulse  response  has  been  shown  to  offer  a  physically  intuitive 
view  of  some  relatively  complex  scattering  phenomena.  The  waveforms 
presented  in  Chapter  IV  are  the  most  accurate  measured  impulse  responses 
available  to  date.  As  a  first  look,  it  is  only  natural  that  they 
present  more  questions  than  answers.  But  the  intimate  relation  between 
the  transient  waveform  and  the  scattering  geometry  directs  further 
research  toward  specific  scattering  mechanisms.  In  keeping  with  the 
straightforward  nature  of  the  impulse  response,  the  scattering  analysis 
has  been  approached  from  an  intuitive  point  of  view. 

We  have  seen  that  the  early  time  behavior  of  the  impulse  response 
waveform  shows  a  clear  and  predictable  connection  to  surface  curvature 
at  the  specular  point.  For  those  aspects  at  which  the  scatterer  has  a 
reflected  field  the  specular  impulse  is  insensitive  to  the  polarization 
of  the  incident  field.  This  is  consistent  with  physical  optics.  But 
the  return  for  times  immediately  following  the  specular  Impulse  is 
polarization  dependent.  The  relation  between  polarization,  surface 
curvatures  in  the  principal  planes,  and  the  impulse  response  wavefonn 


74 


has  been  discussed  for  specific  examples.  For  those  aspects  at  which 
there  is  no  reflected  field  the  initial  response  is  quite  sensitive  to 
polarization. 


The  late  time  behavior  has  displayed  some  characteristic  low 
frequency  resonant  behavior  in  a  variety  of  forms.  The  measured 
waveforms  may  support  efforts  to  extract  natural  poles  by  analytical 
methods  for  the  purpose  of  target  identification  [25],  The  measurements 
show  that  not  only  low  frequency  transients  appear  in  the  late  time,  but 
also  high  frequency  multiple  diffraction  and  creeping  wave  mechanisms. 

Finally,  by  displaying  several  type  of  time  dispersion,  the 
measured  waveforms  suggest  ways  of  extending  the  Uniform  Theory  of 
Diffraction  solution  to  account  for  lower  frequency  phase  dispersion 
associated  with  caustic  regions. 
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APPENDIX  A 


UTD  AND  PHYSICAL  OPTICS  APPROXIMATIONS 

The  physical  optics  approximation  to  the  impulse  response  of  a 
scatterer  is  derived  In  [4]  as, 

F  (t)  -  -1  d2A(Z)  ,  Z  *  ct  .  (A.l) 

1  T 

where  A(Z)  is  the  cross  sectional  area  Intercepted  by  a  plane  at  Z 

parallel  to  the  xy  plane  and  moving  in  the  positive  Z  direction  at  a 

speed  c/2.  The  form  of  A(Z)  and  its  second  derivative  is  shown  In 

Figure  A-l  for  the  sphere.  Figure  A-2  for  the  2:1  right  circular 

cylinder  at  broadside,  and  Figure  A-3  for  the  disc  of  radius  a. 

For  the  broadside  cylinder  In  Figure  A-2  the  delta  function  at  Z*=0 

is  not  well  defined.  While  the  shape  of  this  curve  is  accurate  the 

proper  weight  of  the  impulse  would  require  a  limiting  procedure  for  the 

cylinder  slightly  off  broadside  as  Is  done  in  the  UTD  expression  of 

+  3/2 

Equation  (4.13).  The  behavior  near  Z  ♦  0  varies  as  (*.)  indicating 
an  Ill-defined  Laplace  transform. 

In  Figure  A-3  d%V(Z)  for  the  disc  has  the  form  0f  a  shifted  tangent 
dZ* 

curve.  In  the  limit  as  9  +  u/2  the  curve  approaches  a  doublet. 
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The  first  order  UTO  solution  for  the  sphere-capped  cylinder  Is  well 
documented  In  T.T.  Chu's  dissertation  [12].  A  sketch  of  the  scattering 
geometry  Is  given  In  Figure  A-4.  For  0  <  $  <  90°,  there  is  a  reflected 
field  from  the  sphere,  a  diffraction  from  the  discontinuity  In  surface, 
curvature  at  the  junction  between  sphere  and  cylinder,  and  a  diffraction 
from  the  rear  edge.  For  90  <  $  <  180°,  there  Is  diffraction  from  two 
stationary  points  on  the  end  of  the  cylinder  and  a  diffraction  from  the 
sphere-cylinder  junction. 
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